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We demonstrate, in the context of quadratic fermion lattice models in one and two spatial di- 
mensions, the potential of entanglement renormalization (ER) to define a proper real-space renor- 
malization group transformation. Our results show, for the first time, the validity of the multi-scale 
entanglement renormalization ansatz (MERA) to describe ground states in two dimensions, even at 
a quantum critical point. They also unveil a connection between the performance of ER and the log- 
arithmic violations of the boundary law for entanglement in systems with a one-dimensional Fermi 
surface. ER is recast in the language of creation/annihilation operators and correlation matrices. 
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The renormalization group (RG), concerned with the 
change of physics with the observation scale, is among the 
main ideas underlying the theoretical structure of statis- 
tical mechanics and quantum field theory, and of central 
importance in the modern formulation of critical phe- 
nomena and phase transitions [1]. Its influence extends 
well beyond the conceptual domain: RG transformations 
are also the basis of numerical approaches to the study 
of strongly correlated many-body systems. 

In a lattice model, a real-space RG transformation pro- 
duces a coarse-grained system by first joining the lattice 
sites into blocks and then replacing each block with an ef- 
fective site [2] . Two very natural requirements for a such 
RG transformation are: (i) it should preserve the long- 
distance physics of the system; (ii) when this physics is 
invariant under changes of scale, the system should be a 
fixed point of the RG transformation. 

For the important case of a quantum system at zero 
temperature, the first requirement is fulfilled if, as de- 
termined by White in his density matrix renormaliza- 
tion group (DMRG) [3], the vector space of the effective 
site retains the local support of the ground state. En- 
tanglement renormalization (ER) [4] has recently been 
proposed in order to simulatenously meet the second re- 
quirement. By using disent anglers, ER aims to produce a 
coarse-grained lattice locally identical to the original one, 
in the sense that their sites have the same vector space 
dimension. When this is accomplished, the original sys- 
tem and its coarse-grained version can be meaninglfully 
compared, e.g. through their Hamiltonians or ground 
state properties, leading to a proper real-space RG flow. 

Promisingly, ER has been successfully demonstrated 
for the ID quantum Ising model with transverse mag- 
netic field, where it has been shown that, indeed, at the 
quantum critical point the system is invariant under the 
resulting RG transformation However, plenty of work 
is still required to characterize the main features and 
range of applicability of this new approach and, in par- 
ticular, it remains to be seen whether ER also works in 
the computationally more challenging context of 2D lat- 
tice systems, where DMRG can no longer analyse large 



systems. 

In this paper we explore the performance of ER in 
systems of spinless fermion both on ID and 2D lattices, 
as specified by the quadratic Hamiltonian 

N ^ N 

H = ^ [ala s + 7 (a^al + a s a r )] - A ^ a£a r , (1) 

r,s=l r—1 

where A and 7 are the chemical and pairing potentials and 
the first sum involves only nearest neighbors. In spite 
of its simplicity, Hamiltonian H contains a rich phase 
diagram as a function of A and 7, including insulating, 
conducting and superconducting phases [s[ . Importantly, 
the corresponding ground states span all known forms of 
entropy scaling [a, la]- In addition, H can be diagonalized 
through linear (Fourier and Bogoliubov) transformations 
of the fermion operators a and a* while, by Wick's the- 
orem, all properties of its gaussian ground state |^ GS ) 
can be extracted from the two-point correlators (a\a s ^ 
and (a r a s ). Then, provided that our RG transformation 
also maps fermion modes linearly, the entire analysis can 
be conducted in the space of two-point correlators and 
quadratic Hamiltonians of N fermionic modes, as repre- 
sented by N x N matrices. Hence quadratic fermionic 
models such as (pQ) offer an appealing testing ground for 
ER, one where computational costs have been greatly 
simplified (e.g. H can be diagonalized exactly with just 
0(N 3 ) operations) while keeping a rich variety of non- 
trivial ground state structures. 

We start by rephrasing, in the language of correlation 
matrices, the process of coarse-graining a D-dimensional 
(hypercubic) lattice. We assume that the system is in the 
ground state |^ GS ) of H, which we compute using stan- 
dard analytic techniques (see e.g. [5]). It is convenient 
to redraw the hypercubic lattice so that each site con- 
tains P = p D fermion modes for some integer p. Then a 
hypercube of 2 D sites defines a block that contains P2 D 
modes. The goal of the RG transformation is to replace 
this block with just one effective site made of P' modes, 
with P' < P2 D . We would like to have P' = P, so that 
the sites of the coarse-grained and original lattices are 



identical and we can compare the corresponding Hamil- 
tonians or ground-state reduced density matrices. How- 
ever, in the coarse-graining step only modes of the block 
that are disentangled from the rest of the system can be 
removed (see appendix A). As a result, P' often must be 
larger than P. 

For the sake of simplicity, we continue the analysis for 
the case of a ID lattice (see appendix B for construction 
deatils of MERA for the 2D lattice). Let us temporarily 
replace the N spinless fermion operators a in Eq. (pQ) 
with 27V (self- adjoint) Majorana fermion operators c, 



c 2r -i = a r + al, 



7.t 



C2r 



(2) 



The ground state |^ G s) is then completely specified by 

(c r c s ) = S rs -\- iY rsj (3) 

where T, henceforth referred to as the correlation ma- 
trix, is real and antisymmetric. Similarly, the reduced 
density matrix p GS for a block made of 2 sites, that is 
with L = 2P spinless modes (equivalently, 2L Majorana 
modes) is described by a 2L x 2L submatrix Yl of T. 
This matrix is brought into (block) diagonal form by a 
special orthogonal transformation V, 
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FIG. 1: Top: A block of two sites (four modes) is coarse- 
grained into an effective site by first applying disentanglers U 
across the boundary of the block and then using isometry W 
to project out two modes. Bottom: Same RG transformation 
written in the language of correlation matrices, Eq. J7}. 
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where R G SO(2L) and 



where < v r < 1 are the eigenvalues of Yl, each one 
associated with a pair of Majorana fermions. These pairs 
recombine into L spinless fermions in a product state [9[ 
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where g ri the state of a spinless fermion mode, is mixed 
if v r < 1 and pure ifv r = 1. Notice that since the ground 
state |^gs) is a pure state, a mode in a mixed state must 
be entangled with modes outside the block, whereas a 
mode in a pure state is unentangled from the rest of the 
system. We build an effective site by removing from the 
block, or projecting out from Tl, all the modes that are 
unentangled (pure), and just keeping those P' modes that 
are entangled (mixed). In this way, the coarse-grained 
lattice retains the ground state properties, see appendix. 

The key idea of ER, see Fig. ([1]), is to use disentangling 
unitary transformations, or disentanglers, to diminish P' 
by increasing the number of modes in the block that are 
unentangled from the rest of the system. A disentan- 
gler is implemented through a special orthogonal matrix 
U G SO(2L) that acts on two neigboring sites across the 
boundary of the block, wheareas the coarse-graining is 
implemented by an isometry W = RYp> that selects the 
P' spinless fermion modes to be kept in the effective site, 
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Let T3L describe three consecutive blocks. Then the cor- 
relation matrix T' L for the effective site reads (Fig. ([!])) 



= w^{u@ u) f r 3L (u © co w. 



(7) 



Similarly, the correlation matrix Y' L for the modes to be 
removed is 



T' L = W*(U® U) ] T 3L (U © U) W, 
W = R(Y L -Y P ,), 
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Our goal is to maximize the purity of the modes to 
be projected out, so that they become as unentangled as 
possible. The sum of their purities, Y^r=P'+i v n 18 nan ° °^ 
the antisymmetric trace of T' L , tr(f^Y^). Consequently, 
U and W are obtained from the optimization 

tr (f' L Y[) , (10) 



max 

U,ReSO(2N) 



that we address through a sequence of alternating opti- 
mizations for U and R [111 ]. 

Then, given the correlation matrix T for |\I/gs)? the 
RG transformation is implemented in three steps: (i) 
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FIG. 2: Scaling of the entanglement entropy Sl 9] in ID sys- 
tems. Left: Quantum Ising model, 7=1. Bold (solid/dotted) 
lines represent entanglement at criticality, A = 1. The sys- 
tem is an entangled fixed point of our RG transformation: 
the correlation matrices {T^, T^ 2 \ • • • } quickly converge to a 
fixed r* sing . In particular, the renormalized entanglement of 
a block is constant. Thin lines correspond to a non-critical 
system, A = 1.001, which the RG flow eventually brings a 
product (unentangled) ground state. Right: Quantum XX 
model, 7 = 0. Bold/thin lines represent two critical cases, 
A = and A = cos(157r/16). They belong to the same uni- 
versality class and are found to indeed converge to the same 
correlation matrix Txx 5 (with Tx X 7^ r* sing ) and in particular 
to the same renormalized entropy. 



first a submatrix Tsl for three consecutive blocks is ex- 
tracted from T; (ii) then disentangler U and isometry 
W are computed using the optimization ([TP]) while keep- 
ing P' = P modes in the effective site; (Hi) finally, U 
and W are used to transform the original TV-mode sys- 
tem into a coarse-grained system with just N/2 modes 
and effective correlation matrix I^ 1 ). Some of the modes 
that are removed are still slightly mixed. Their mixness 
e r = 1 — v r quantifies the errors introduced. Iteration of 
the RG transformation produces a sequence of increas- 
ingly coarse-grained lattices, described by correlation 
matrices {T^ 1 ), T^ 2 \ • • • }. The corresponding disentan- 
glers {U^\U^ 2 \---} and isometries {W^\ W^ 2 \ ■ ■ ■ } 
constitute the multi-scale entanglement renormalization 
ansatz (MERA) |7[ for the ground state |^ G s)- 

We have applied the present RG approach to Hamil- 
tonian ([T]) in the thermodinamic limit N — > oo. First 
we consider ID systems, where the whole (7, A) plane 
can be mapped into the quantum spin XY model us- 
ing a Jordan- Wigner transformation, (i) For the line 
7 = 1 [equivalent to the quantum spin Ising model] we 
consider P = 2 modes per site and apply 13 iterations 
of the RG transformation, so that a final effective site 
(with just P = 2 modes) corresponds to 2 x 2 13 = 16384 
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FIG. 3: Dispersion relation of Hamiltonian JTJ) in ID with 7 = 
0, quantum spin XX model, under successive RG transforma- 
tions. Shading indicates the Fermi sea. A sequence of lo- 
cal, coarse-grained Hamiltonians is obtained {H^\ H^ 2 \ • • • } 
with their corresponding dispersion relations {z/i, z/2, • • • } con- 
verging to a straight line, a fixed point of the RG flow. Con- 
vergence is achieved very quickly at half filling (A = 0) and 
slower for A = cos(37r/4). These results have been obtained 
by minimizing the energy (Sect. IV of Ref. [l(j) while keeping 
8 modes in each effective site. 



modes of the original system. At the critical point A = 1, 
which is the most demaning, the mixness of the removed 
modes is at most e r = 1.2 x 10 -4 . The effect on local 
observables, even after the 13 iterations, is remarkably 
small: the error in the critical ground state energy is 
less than 10 _T , while the two-point correlators (a£a s ), 
reconstructed from the MERA, accumulate a relative er- 
ror that ranges from 10 _T for nearest neighbours to 10% 
for \r — s\ w 4,000. Had we not used disent anglers, the 
error in the energy would be 10 -3 after only a single RG 
transformation and an error of 10% in the two-point cor- 
relators is already achieved for \r — s\ = 42. (ii) The 
line 7 = [equivalent to the quantum spin XX model] is 
critical for |A| < 1. Here we consider P = 4 modes per 
site and apply again 13 iterations of the RG transforma- 
tion, reaching sizes of 4 x 2 13 = 32768 modes. The errors 
in energy and correlators are similar to those in the line 
7 = 1. In both cases, an analysis of the RG flow and 
its fixed points in terms of entanglement is quite insight- 
ful, see Fig. ([2]). ER can also be used to generate a RG 
transformation in the space of Hamiltonians, by replac- 
ing Eq. ([TO]) with a minimization of the energy (see Sect. 
IV of Ref. for details). Fig. ([3]) shows that critical 
systems are also fixed points of this alternative approach, 
that preserves the low energy spectrum. 

In 2D the model has three phases, denoted I, II and 
III in Ref. where the distinct forms of entanglement 
scaling were characterized. In phases II (critical, with a 
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Block length, L 

FIG. 4: Entanglement entropy .Sl of a block of L x L modes in 
2D models. Left: in the critical phase II and the non-critical 
phase III (bold/fine lines respectively) the entanglement en- 
tropy grows linearly with the size L of the boundary of the 
block, Sl ~ L (boundary law). As in ID, the renormalized 
entanglement is constant for the critical model and it eventu- 
ally vanishes for the non-critical model. We have considered 
7 = 1 and A = 2, A = 2.05 for the critical/non-critical case. 
Right: The critical phase II system (7, A) = (1,2) is replotted 
for comparison against critical phase I, (7, A) = (0,0), where 
the system has a ID Fermi surface and the entanglement en- 
tropy has a logarithmic correction, Sl 00 LlogL. Here disen- 
tanglers are not able to reduce the renormalized entanglement 
down to a constant. 



Fermi surface consisting of a finite number of points) and 
III (non-critical, with a gap in the energy spectrum) we 
are once more able to coarse-grain the system in a quasi- 
exact, sustainable manner. This is remarkable. The en- 
tropy of a square block made of L 2 modes grows as the 
size of its boundary, Sl ~ L This implies that the 
number of modes we should keep in an effective site grows 
exponentially with the number of iterations of the RG 
transformation, which is precisely why DMRG does not 
work for large 2D systems. Instead, disentanglers bring 
this number again down to just a constant. As a re- 
sult one can, in principle, explore systems of arbitrary 
sizes. In particular, by considering P = 4 2 modes per 
site we apply r = 4 iterations of the RG transformation, 
with a final block effectively spanning P x 4 r+1 = 16384 
modes, whilst maintaining truncation errors of the same 
scale as the ID models analysed, e r = 1.1 x 10 -4 . As 
in the ID case, the structure of fixed points of the RG 
flow can be understood in terms of the renormalized en- 
tanglement, see Fig. (j4]). On the other hand, Phase I 
(critical, with a one-dimensional Fermi surface) is so en- 
tangled that ER is no longer able to prevent the growth 
in the number P' of modes that need to be kept per site. 
The system displays a logarithmic correction to the en- 
tropy, S l ~ £ log L [1, @, H, while the MERA can only 



reproduce a linear scaling Sl ~ L [7] if just a constant 
number of modes are kept per site, P' = P. 

We have presented, in the simplified context of fermion 
models with quadratic Hamiltonian, unambiguous evi- 
dence of the validity of the ER approach in ID and 2D 
systems. Similar derivations can be also conducted for 
bosonic lattice systems with quadratic Hamiltonians [l2[ • 
Our results show, for the first time, that the MERA 
Q is an efficient description of 2D ground states. A 
number of examples also confirm that (i) ER produces 
a quasi-exact, real-space RG transformation where the 
coarse-grained lattice is locally equivalent to the original 
one, enabling the study of RG flow both in the space of 
ground states and Hamiltonians; (ii) non-critical systems 
end up in a stable fixed point of this RG flow, where the 
corresponding ground state is a product (i.e. fully dis- 
entangled) state, whereas critical systems end up in an 
unstable fixed point, with an entangled ground state. 

Interestingly, ER also sheds new light into the ground 
state structure of two-dimensional systems with a one- 
dimensional Fermi surface: only when logarithmic correc- 
tions appear for the entropy Sl of a large block [H|, 0, Hj] , 
does the above simple picture break down, suggesting the 
need for a generalized MERA [12| in order to describe 
such systems. 
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Appendix A. — Here we describe the process of 
coarse-graining a lattice by replacing blocks of sites with 
effective sites. We show that the truncation of the Hilbert 
space of a given block can be implemented by eliminating 
some of the modes in that block. 

We consider a fermionic lattice system in its (gaussian) 
ground state |\£ G s)- Let V be the vector space of a block 
containing L modes and let p GS denote the reduced den- 
sity matrix of |^) GS on the block. We assume that the 
support of \SP) GS is concentrated in a subpace V p C V. 
Then, following White 0, the optimal coarse-graining of 
the block is obtained by defining an effective site s' with 
vector space V s = V p . In our case, p GS is the tensor 
product of density matrices g r for individual modes 

L L / 1+V r n \ 

Pos = ®0r = ®( I ( n ) 

Suppose that the first P' modes are in a mixed state and 
the remaining L — P' modes are in a pure state. Then 
we can write 

p' L 

Pgs = ((g) Qr) ® ( (g) Q r )=(J^lT, (12) 
r=l r=P'+l 

where a is a mixed state with rank 2 P whereas tt is 
a projector with rank 1. Let V = Y a V 7 *" be a ten- 
sor factorization of V such that a = try^(p). The key 
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observation is that Y a = V p , and that p and a have the 
same none- vanishing eigenvalues. Therefore, we have two 
equivalent ways of constructing the space V s for the ef- 
fective site s' while preserving the support of the ground 
state density matrix p GS . On the one hand, V s can be 
obtained by projecting V on the support V p of p. On 
the other, V can also be build by factorizing the space V 
into two factor spaces Y a and V 71 ", and by then tracing 
out the second factor, corresponding to modes in a pure 
state. Both constructions lead to an equivalent effective 
lattice. Finally, tracing out the factor space of mode r 
corresponds, in the language of correlation matrices T^, 
to removing the rth row and a rth column of VTlV^ in 
Eq. process to which we referred to as projecting out 
the mode. 

Appendix B. — Here we describe the multi-scale 
entanglement renormalization ansatz (MERA) that we 
have used in the present work for 2D systems. It differs 
from the one described in Refs. 0,0]. 

The MERA can be understood as a peculiar class of 
quantum circuit with bounded size causal cones 0, 0- 
The causal cone structure of the circuit is the key prop- 
erty that allows for efficient computations with this 
ansatz. In practical realizations, the detailed local struc- 
ture of the MERA (how disentanglers and isometries are 
interconnected) depends on the specific problem under 
consideration. In particular, for a 2D system the geom- 
etry of the lattice (square, triangular, Kagome,...) or 
the symmetries of the state that we intend to represent 
will be taken into consideration in order to choose a spe- 
cific realization of the MERA. Fig. ([5j) describes the 2D 
MERA for a square lattice discussed in [4|, [7[ . Disentan- 
glers and isometries act and reduce one lattice direction 
(say x or y) at a time. This seems to be the most eco- 
nomical realization of a 2D MERA in terms of how the 
computation cost scales with the index dimension \. 

Fig. ([6]) describes instead the realization of the 2D 
MERA used in the present work. In this case both lat- 
tice directions are addressed simultaneously. While this 
realization of the MERA induces a cost that scales as a 
larger power of x, it seems to be more adequate for prob- 
lems where it is important to preserve symmetry under 
90° rotations. Similarly, one could consider other realiza- 
tions, provided the fundamental causal cone properties of 
the associated quantum circuit are preserved. 

NOTE ADDED: The 2D MERA realization used in 
this work has also been subsequently but independently 
discussed in Ref. [13]. Appendix B has been added to 
clarify Jjow our present approach differs from the one in 
Refs. 
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FIG. 5: 2D MERA on a square lattice as described in j3,Q. 
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FIG. 6: 2D MERA on a square lattice as used in the present 
calculations. 
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